Motion of Vacancies in a Pinned Vortex Lattice: Origin of the 

Hall Anomaly 

P. Ao 

Department of Theoretical Physics, Umea University 
00 ■ 901 81, Umed, Sweden 

0\ . (March 3, 2009) 



O 

(N 



a 



1 



Abstract 



^ . Physical arguments are presented to show that the Hall anomaly is an effect 

OO _ 
^^ . of the vortex many-body correlation rather than that of an individual vor- 

(N 

O 

0> , vacancy formation energy, is obtained for thin films. At low temperatures a 



tex. Quantitatively, the characteristic energy scale in the problem, the vortex 



scaling relation between the Hall and longitudinal resistivities is found, with 
the power depending on sample details. Near the superconducting transition 



c . 

O . temperature and for small magnetic fields the Hall conductivity is found to 

be proportional to the inverse of the magnetic field and to the quadratic of 
the difference between the measured and the transition temperatures. 
PACS#s: 74.60.Ge 
The ubiquitous occurrence of the Hall anomaly in the mixed state of both conventional 
and oxide superconductors, the sign change of the Hall resistivity below the superconducting 
transition temperature and the smallness of the Hall angle, has defied a consistent expla- 
nation so far 1^. A straightforward application of the Magnus force cannot explain this 
phenomenon. This failure leads to a frustration of questioning the Magnus force as the only 
transverse force in the vortex dynamic equation. The transverse force has been subsequently 
modified into various forms [Q]. Following the general properties of a superconductor, on 
the other hand, recent theoretical work on vortex dynamics have shown that there is no 
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other transverse force than the Magnus force, a result of the topological constraint and the 
momentum conservation. |^Q One can further show that the Magnus force is equivalent to 
the spectral flow process. An apparent conflict between the theoretical reasonings and 
the experimental measurements naturally arises. 

In the present letter we attempt to solve this puzzle by showing that the Hall anomaly 
can be understood based on the vortex vacancy motion in a pinned vortex lattice, and 
emphasize that the anomaly is a property of the vortex many-body correlation rather than 
that of an individual vortex. We will demonstrate that vacancies can have the lowest energy 
scale, and they dominate the thermal activation contributions to the vortex motion at low 
temperatures. The present vacancy model for the Hall anomaly is also consistent with other 
measurement such as the Nernst effect. In the following we present our arguments leading to 
the model, and discuss its predictions and validity conditions. For simplicity, we will consider 
an isotropic s-pairing superconductor with one type of charge carriers in two dimension. In 
this situation vortices(or straight vortex lines) can be viewed as point particles. |^ 

The vortex dynamic equation for a j-th vortex of unit length in the sample takes the 
form of the Langevin equation identical to that of a charged particle in the presence of a 
magnetic fleld: 

mi-j = q-^h (v^,j - Tj) x z - titj + Fp + f , (1) 

with an effective unit length mass m, a pinning force Fp, a vortex viscosity rj, and a fluc- 
tuating force f. In Eq.(l) q = ±1 is the vorticity, h the Planck constant, ps the superfluid 
electron number density at temperature T, and z the unit vector in 2;-direction. If there is 
a temperature gradient, the thermal force F^ = —S(f,VT should be added in at the right 
hand of Eq.(l), with s^ the entropy carried by a vortex. The term associated with the total 
superconducting electron velocity v^ t = v^ + Vs,in and the vortex velocity r at the right side 
of Eq.(l) is the Magnus force. Although v^^j is due to all other vortices, here we split it 
into two parts, with v^ corresponding to the rearrangement of vortices due to the externally 
applied supercurrent and Vs,in accounting for the rest contribution describing the vortex 



interaction without external current. In the following we will assume v^ is small such that 
this splitting is valid. 

It is evident that in the mixed state of any real superconductor the many-body correlation 
between vortices and the pinning effect cannot be ignored. The competition between them 
is the source of the rich static and dynamical properties of flux phases []^,|^ . We will take the 
Abrikosov lattice as the known manifestation of the many-body correlation for the starting 
point to advance our arguments. The vortex pinning is also important in our reasoning, 
though several quantitative results obtained below do not explicitly depend on it. If there 
were no pinning for vortices, the whole vortex lattice would move together under the influence 
of an externally applied current in the same manner as that of independent vortices. Hence 
one would get the same sign of the Hall resistivity in both superconducting and normal 
states. In the presence of pinning centers in the sample as well as the edge pinning, the 
vortex lattice will be pinned down. In such a situation the motion of the vortex lattice is 
made possible by various kinds of defect motions due to thermal fluctuations. We will argue 
below that at low temperatures the dominant contribution to the motion is due to vortex 
vacancies, and the Hall anomaly occurs. 

For two vortices separated by a distance r, which is less than the effective magnetic 
screening length A^ = X\/d{d < \l, X± = \l H d > Xl) but greater than ^o, the interaction 
potential is Vo{r) = 2d($o/47rAL)^ ln(r/^o)- Here A| = m*c^/Mpse^ is the London 
penetration depth, m* the effective mass of a Cooper pair, ^o the coherence length of the 
superconductor, and d the thickness of the superconductor film. The energy scale eo = 
(i($o/4vrAL)^ sets both the scale for the strength of vortex interaction and the scale for the 
strength of a strong pinning center. The energy for a dislocation pair separated by a distance 
larger than the vortex lattice constant is given by Vd{r) = (eo/2v^7r) ln(r/ao), ^,0 with Oq 
an order of the the vortex lattice constant. The energy scale eo/2v^vr for the dislocation 
pair here is about 10 times smaller than eo for the vortex interaction and pinning centers. 
It is energetically favorable to have dislocation pairs in the lattice. Hence for temperature 
T << eo we can ignore the contribution from the vortices hopping out of pinning centers and 



the creation of vortex-antivortex pairs. The vortex lattice is then effectively pinned down. 
Because vacancies and interstitials can be viewed as the smallest dislocation pairs |T0|, we 
immediately have the estimated energy scale for vacancy formation energy e^ as, by putting 
r ~ 2ao in Vd{r), 

1 ^*«^^.. (2, 
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This result is valid for an intermediate magnetic field B: Hd < B < Hc2/'2. For thicker films 
the thickness in Eq.(2) will be replaced by a crossover thickness dc due to the 2;-direction 
correlation, whose precise value is a complex and unknown function of various parameters 
such as the magnetic field, the pinning, the temperature, and anisotropy. In the case dc 
is finite, its estimation in the high magnetic field limit is as follows. Ring type vacancy 
excitations are possible in thicker films. Its energy scale is determined by the smallest 
ring, which should be the size of the vortex lattice constant. In this case dc is an order 
of the vortex lattice constant, dc ~ clq- However, because of the large anisotropy in the 
HTcS materials, dc can be the order of the CuO layer spacing close to the superconducting 
transition temperature T^q. 

It is clear from the above analysis that vacancies and interstitials have the lowest ex- 
citation energy scales. We note that the value in Eq.(2) is consistent with the variational 



and numerical calculations 111, with about a factor 2 smaller, and also with the estimation 



from the dislocation core energy [|T^. Using the shear modulus results for B < Hd ^ and 
Hc2/2 < B < Hc2 00 we have obtained the corresponding vacancy formation energies as 
(2/97r)i/2(6o/AL)^/2e"''o/^^ eo with b^ = %/B, and 0.7/V37r(l - B/Hc2f eo, respectively. 

Now we argue that the vacancy formation energy is even lower than that of an inter- 
stitial. The experimental observations at low magnetic fields have shown the abundance 



of vacancies comparing with interstitials. |T3[ The natural explanation is that the the va- 



cancy formation energy is lower than that of interstitials, therefore by thermal fluctuations 
vacancies have a higher density. The theoretical calculations have also confirmed the lower 
vacancy formation energy [jl4|. This phenomenon of the vacancy formation energy is lower 



than that of interstitials has also been observed in other crystalhne structures [0. We 
conclude that vacancies will dominate thermal fluctuation contributions to resistivities at 



low enough temperatures for low magnetic fields. |T^ 

We show next that in the pinned vortex lattice a vacancy behaves as a vortex with a 
vorticity —q and an interstitial as a vortex with +g, respectively. Let u be the displacement 
vector at position r, with a point defect, vacancy or interstitial, at fq. According to Eq.(l) 
the transverse force acting on the defect is, measured from the pinned perfect vortex lattice. 



Tpd 



qPlh Jdh6p{r) (v,-u)xi. (3) 



Here the vortex density 6p deviated from a perfect lattice is determined by the dilatation 
V ■ u: 5p = V • u/Sq, with 5*0 the area of a unit cell in the vortex lattice. By definition, 

V-u = TSoS\r-ro), (4) 

with ' — ' for a missing vortex, a vacancy, and '+' for an extra vortex, an interstitial. Using 
Eqs.(3,4), we have the desired transverse force on the defect as 

Kf = Tq^h (v, - i-o) X z . (5) 

This is identical to the dynamics of a hole or a particle in a semiconductor in the presence 
of a magnetic field, with a pinned perfect vortex lattice as a filled valence band and a 
vacancy in real space as a hole in the energy space. Eq.(5) shows that both a vacancy and 
an interstitial will move along the direction of the applied supercurrent v^. This implies 
that vortices defining vacancies move against the direction of v.,, a result of the many- 
body correlation and pinning. This leads us to our main conclusion that at low enough 
temperatures the sign of the Hall resistivity is different from its sign in the normal state 
because of the dominance of vacancies. Quantitatively, vacancies and interstitials may be 
considered as independent particles moving in the periodic potential formed by the vortex 
lattice and a random potential due to the residue effect of pinnings. The potential height of 
the periodic potential as well as that of the random potential is an order of e„. Assuming 
the vacancy (interstitial) density n^{n.;) in a steady state, the longitudinal resistivity is 



_ _(?_ Y^ Vi psh/2 nj_ 
P^^ - 2e- t, ril + iPsh/2y p. ' ^'^ 

and the Hall resistivity 

_ _/?_ ^ {psh/2f TH 

with Qv = —q and qi = q. Here rj^^i are the effective vacancy and interstitial viscosities, 
related to their diffusion constants in the periodic potential due to the vortex lattice by 
the Einstein relation between the diffusion constant and the mobility. It should be pointed 
out that contributions of other vortex motions to resistivities such as vortex-antivortex 
pairs, which are omitted here for their smaller activation probabilities, are additive to those 
of vacancies, and that the including of the normal fluid (quasiparticle) contributions is 
straightforward [0. 



Under the driving of a temperature gradient, the effective thermal force felt by a vacancy 
is opposite in sign to the force felt by an interstitial or a vortex in direction but equal in 
magnitude, F^ = H-s^VT. This can be seen by repeating the demonstration from Eq.(3) 
to (5). Then the Nernst effect due to vacancies has the same sign as that of vortices or 
interstitials. Therefore our model gives that in the Hall anomaly regime there is no sign 
change for the Nernst effect, and furthermore, the Nernst effect is more pronounced because 
of the additive contributions due to both vacancies and interstitials. This is in agreement 
with the experimental observations [|l^]. 

Before exploring of consequences of Eqs.(6,7) we discuss the qualitative implications of 
the present model. In the above picture, to obtain a maximum contribution of vacancies, 
we need the vortex lattice to define vacancies and a sufficiently strong pinnings to prevent 
the sliding of vortex lattice. The existence of a whole lattice structure is nevertheless un- 
necessary. Sufficiently large local crystalline structures, like lattice domains, will be enough 
to define vacancies. Therefore vacancy-like excitations in a vortex liquid state can exist, 
because of the presence of large local orderings. Whether or not this is also true for a vor- 
tex glass state depending on details. For example, a further lowering of temperature may 



quench a vortex system into a glass state with no local crystalline structure. Then vacancies 
will disappear and the sign of the Hall resistivity will change again. On the other hand, for 
a fixed temperature if the pinning is too strong, for example, the (random) pinning center 
density is much larger than the vortex density, vortices will be individually pinned down 
and the local lattice structure required for the formations of vacancies and interstitials will 
be lost. This suggests that the Hall anomaly only exists in a suitable range of pinnings 
and magnetic fields, that is, for Bi < \B\ < Bu with the lower and upper critical fields 
determined by pinning as well as by temperature. 

Now we study the limiting cases of Eqs.(6,7). At low temperatures the motions of 
vacancies and interstitials in the vortex lattice are thermal hopping: r]^ = rjQ e"" '^^I^bT ^t^^ 
Vi = Vo e"' '^■"/^bT^ with a„, aj(presumably a„ < ai) numerical factors of order unity and r]Q 
insensitive to temperature. In this limit, the vacancy (interstitial) density n^, = uq e^^" ^v/ksT 
{rii = no e~^' <^v/kBT^^ with b^ = l(6.j > 1) for the thermally activated vacancies(interstitials) 
and bv{bi) = for the pinning center induced vacancies(interstitials). In the following we 
further assume that i]^,rii » Psh/2, corresponding to the Hall angle | tan^| = \pyx/Pxx\ << 
1 common in experiments. Under this assumption, we obtain the Hall angle as 



tan^ = —q 



(8) 






Here 7 = 2ai + bi — 2a„ — b^. The high temperature limit ksT > e^ is achieved near 
superconducting transition temperature T^o, but the thermal creation of a vortex-antivortex 
pair is still improbable, because the relevant energy scale eo is about 10 times bigger than e„. 
In the low temperature limit both longitudinal and Hall resistivities vanish exponentially. 
We obtain a scaling relation between them as 

Pyx = ^ pL , (9) 

with A = -g(ps/i/2r7o)^"/^""+^''^(2e^Ps//;,no)"''/^""+^''\ and the power 



^ = — ^— ' (10 

varying between 1 and 2, depending on the detail of a sample which determines the numerical 
factors a,^ and h,^. If all vacancies are produced by pinnings, we have 6^, = and v = 2. 
In this case A is independent of B because tiq is. In the other limit, if all vacancies are 
produced by thermal activations, and if a„ << 1, we have 6„ = 1 and z/ ~ 1. In this case A 
will be independent of B if if rj^ is. 

Another useful quantity is the Hall conductivity a^y = Pyx/iplx + Pyx)- Under the 
same assumption oi rj^^rji » Psh/2 we obtain the Hall conductivity due to vacancies and 
interstitials, from Eqs.(6,7), as 

2g2 p g-(2at,+6„)et,/fcsT _ g-(2ai+6,)ei,/A:sT 



<yxy = -q- 



h Hq \^-{av+K)<Sv/kBT _|_ ^-{ai+bi)ey/kBTy 



^ h no A ksT T a — V 

As discussed above, here < 6^, < 1 and 7 ~ 0(1). Near the superconducting transition 
temperature Tco, Ps = Pso(l — T/Tco) and e^ = €^{1 — T/Tco) because of the London 
penetration depth in Eq.(2). We may further assume no = -B/$o, with $0 the flux quantum. 
From Eq.(ll) we obtain 

(^xy = «! ^ , (12) 

with tti = -q{2e^/h)pso%ieyo/4:kBTco. Taking p^o = 10^'^/cm^, 7 = 1, and e^o/^B^co = 50, 
we find |q;i| ~ 20 TfiQ~^cm^^. 

Two comments are in order. 1. A naive accounting of the many-body correlation and 
pinning may not lead to the sign change in the Hall resistivity: Since vortex interaction 
terms cancel each other when summing over all vortices, one would like to conclude that 
there is no many-body correlation effect on the sign of the Hall resistivity. If this claim 
were correct, the same argument would lead to no sign change for the Hall resistivity in a 
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hole semiconductor, and no such phenomena as the quantum Hall effect. This absence of 
Hall anomaly is the result of the underestimation of the many-body correlation. 2. There 
might be a tendency to mix up antivortices and vacancies. As discussed above, the creation 
energy of an antivortex is about 10 times larger than that of a vacancy, which makes it 
energetically unfavorable. Furthermore, since an antivortex feels the same thermal force as 
a vortex, F^ = —s^VT, it has an opposite sign contribution to that of a vortex for the 
Nernst effect, in conflicting with experiments. ^ 

In conclusion, we have demonstrated that within the vortex dynamics equation the Hall 
anomaly can be explained. What has been missed in previous models is a proper considera- 
tion of the competition between the many-body correlation and pinning. We have proposed 
the model of vacancy motion in a pinned lattice as a concrete realization: the characteristic 
energy in the model, the vacancy formation energy, is obtained; and vacancies move along an 
applied supercurrent as the origin for the Hall anomaly. The model leads to an exponential 
tail and the scaling relation at low temperatures, and no sign change for the Nernst effect. 
Near the superconducting transition temperature and for small magnetic fields the Hall con- 
ductivity is found to be proportional to the inverse of the magnetic field and is quadratic in 
the temperature different from the transition. For thin enough films the activation energy 
in the low temperature limit has a linear film thickness dependence. 
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